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We describe the dynamical behavior of both class A and class B lasers in terms of a Lyapunov 
potential. For class A lasers we use the potential to analyze both deterministic and stochastic 
dynamics. In the stochastic case it is found that the phase of the electric field drifts with time in 
the steady state. 

For class B lasers, the potential obtained is valid in the absence of noise. In this case, a general 
expression relating the period of the relaxation oscillations to the potential is found. We have 
included in this expression the terms corresponding to the gain saturation and the mean value of the 
spontaneously emitted power, which were not considered previously. The validity of this expression 
is also discussed and a semi-empirical relation giving the period of the relaxation oscillations far 
from the stationary state is proposed and checked against numerical simulations. 
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I. INTRODUCTION 

Even for non-mechanical systems, it is occasionally 
possible to construct a function (called Lyapunov func- 
tion or Lyapunov potential) that decreases along trajec- 
tories jlj. The usefulness of Lyapunov functions lies on 
the fact that they allow an easy determination of the 
fixed points of a dynamical (deterministic) system as the 
extrema of the Lyapunov function as well as determin- 
ing the stability of those fixed points. In some cases, the 
existence of a Lyapunov potential allows an intuitive un- 
derstanding of the transient and stationary trajectories 
as movements of test particles in the potential landscape. 
In the case of non-deterministic dynamics, i.e. in the 
presence of noise terms, and under some general condi- 
tions, the stationary probability distribution can also be 
governed by the Lyapunov potential and averages can be 
performed with respect to a known probability density 
function. The aim of this work is to construct Lyapunov 
potentials for some laser systems. We start, then, by 
briefly reviewing the main features of the laser as a dy- 
namical system. 

A laser has three basic ingredients: i) a gain medium 
capable of amplifying the electromagnetic radiation prop- 
agating inside the cavity, ii) an optical cavity that pro- 
vides the necessary feedback, and iii) a pumping mech- 
anism. A complete understanding of laser dynam- 
ics is based on a fully quantum-mechanical description 
of matter-radiation interaction within the laser cavity. 
However, the laser is a system where the number of 
photons is much larger than one, thus allowing a semi- 
classical treatment of the electromagnetic field inside 
the cavity through the Maxwell equations. This fact 
was introduced in the semiclassical laser theory, devel- 
oped by Lamb pLpf and independently by Haken [pHTJ. 



This model for laser dynamics was constructed from the 
Maxwell-Bloch equations for a single-mode field interact- 
ing with a two-level medium. The semiclassical laser the- 
ory ignores the quantum-mechanical nature of the elec- 
tromagnetic field, and the amplifying medium is modeled 
quantum-mechanically as a collection of two-level atoms 
through the Bloch equations. A simpler description can 
be obtained by deriving rate equations for the temporal 
change of the electric field (or photons number) inside 
the cavity and the population inversion (carriers number 
in the case of semiconductor lasers) ||. Rate equations, 
with stochastic terms accounting for spontaneous emis- 
sion noise, have been extensively used for semiconductor 
lasers. 

Different types of lasers can be classified according to 
the decay rate of the photons, carriers and material po- 
larization. Arecchi et al. || were the first to use a clas- 
sification scheme: class C lasers have all the decay rates 
of the same order, and therefore a set of three nonlin- 
ear differential equations is required for a satisfactory 
description of the electric field, the population inversion 
and the material polarization. For class B lasers, the po- 
larization decays towards the steady state much faster 
than the other two variables, and it can be adiabatically 
eliminated. Class B lasers, of which semiconductor lasers 
JirH are an example, are then described by just two rate 
equations for the atomic population inversion (or carriers 
number) and the electric field. Another examples of class 
B lasers are CO2 lasers and solid state lasers [jllj . Finally, 
in class A lasers population inversion and material polar- 
ization decay much faster than the electric field. Both 
material variables can be adiabatically eliminated, and 
the equation for the electric field is enough to describe 
the dynamical evolution of the system. Some properties 
of class A lasers, like a dye laser, are studied in 0,[L3|. 
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In this paper we interpret the dynamics of both class A 
and class B lasers by using a Lyapunov potential. 

The paper is organized as follows. In Section || we 
present a brief review of the relation of Lyapunov poten- 
tials to the dynamical equations and the splitting of those 
into conservative and dissipative parts. We consider the 
example of class A lasers. In this case, the Lyapunov 
potential gives an intuitive understanding of the dynam- 
ics observed in the numerical simulations. In the pres- 
ence of noise, the probability density function obtained 
from the potential allows the calculation of stationary 
mean values of interest as, for example, the mean value 
of the number of photons. We will show that the mean 
value of the phase of the electric field in the steady state 
varies linearly with time only when noise is present, in 
a phenomenon reminiscent of the noise-sustained flows. 
In Section III, the dynamics of rate equations for class 
B lasers is presented in terms of the intensity and the 
carriers number (we will restrict ourselves to the semi- 
conductor laser). In this case we have found a potential 
which helps to analyze the corresponding dynamics in the 
absence of noise. By using the conservative part of the 
equations, one can obtain an expression for the period of 
the oscillations in the transient regime following the laser 
switch-on. This expression extends the one obtained in 
a simpler case by an identification of the laser dynam- 
ics with a Toda oscillator in [[uj. Here, we have added 
in the expression for the period the corresponding mod- 
ifications for the gain saturation term and spontaneous 
emission noise. Finally, in section IV, we summarize the 
main results obtained. 



II. POTENTIALS AND LYAPUNOV 
FUNCTIONS: CLASS A LASERS 

The evolution of a system (dynamical flow) can be clas- 
sified into different categories according to the relation 
of the Lyapunov potential to the actual equations of mo- 
tion Jl^,[l6| • We first consider a deterministic dynamical 
flow in which the real variables (xi,...,a;jv) = x satisfy 
the general evolution equations: 



^ = /i(x) ' 



1, 



,N 



(1) 



In the so-called potential Bow, there exists a non- 
constant function V(x) (the potential) in terms of which 
the above equations can be written as: 



dxi 
~~dt 



N 



dV 

j=i ■> 



(2) 



where S(x) is a symmetric and positive definite matrix, 
and Vi (x) satisfy the orthogonality condition: 



N 



dV 

dxi 



(3) 



A non-potential Bow, on the other hand, is one for which 
the splitting (0), satisfying (|J), admits only the trivial 
solution V(x) = constant, Vi(x) = /i(x). 

Since the above (sufficient) conditions for a potential 
flow lead to dV/dt < 0, one concludes that V(x) (when it 
satisfies the additional condition of being bounded from 
below) is a Lyapunov potential for the dynamical system. 
In this case, one can get an intuitive understanding of the 
dynamics: the fixed points arc given by the extrema of 
V(x) and the trajectories relax asymptotically towards 
the surface of minima of V(x). This decay is produced 
by the only effect of the terms containing the matrix S in 
Eq. (||), since the dynamics induced by u, conserves the 
potential, and i>i(x) represents the residual dynamics on 
this minima surface. A particular case of potential flow is 
given when Ui(x) can also be derived from the potential, 
namely: 



A' 



dxj _ ^ n ®V 



j'=i 



dx j 



(4) 



where the matrix -D(x) = 5(x) +A(x), splits into a posi- 
tive definite symmetric matrix, S, and an antisymmetric 
one, A. In this case, the residual dynamics also ceases 
after the surface of minima of V(x) has been reached. 

We now describe the effect of noise on the dynamics of 
the above systems. The stochastic equations (considered 
in the Ito sense) are: 



dxi 
~~dt 



N 



= /i(x)+X)flw(x)6-(t) 

3=1 



(5) 



where <7y (x) are given functions and £j (t) are white noise: 
Gaussian random processes of zero mean and correla- 
tions: 



(6) 



e is the intensity of the noise. 

In the presence of noise terms, it is not adequate 
to talk about fixed points of the dynamics, but rather 
consider instead the maxima of the probability density 
function P(x, t), which satisfies the multivariate Fokker- 
Planck equation JT^JT^] whose general solution is un- 
known. When the deterministic part of ([5]) is a potential 
flow, however, a closed form for the stationary distribu- 
tion P s t(x) can be given in terms of the potential V(x) 
if the following (sufficient) conditions are satisfied: 

1. The fluctuation-dissipation condition, relating the 
symmetric matrix S to the noise matrix g: 



Si 



N 

fe=l 



S = g 



(7) 
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2. Sij satisfies: 



A* 



3'=1 



Vz 



(8) 



This condition is satisfied, for instance, for a con- 
stant matrix S. 

3. u, is divergence free: 



A' 



Edvi 



i=l 



(9) 



this third condition is automatically satisfied for poten- 
tial flows of the form (ji|) with a constant matrix A. 

Under those circumstances, the stationary probability 
density function is: 



P st (x) =Z- l exp[- 



V(x) 



(10) 



where Z is a normalization constant. Graham J19| has 
shown that if conditions || and || arc not satisfied, then 
the above expression for P s t(x) is still valid in the limit 
e -► 0. 

As an example of the use of Lyapunov potentials in a 
dynamical system, we consider class A lasers || whose 
dynamics can be described in terms of the slowly varying 
complex amplitude E of the electric field: 



E=(l 



l + (3\E\ : 



« E + ((t) 



(11) 



where a, f3, T and k are real parameters, k is the cav- 
ity decay rate; T the gain parameter; /3 the saturation- 
intensity parameter and a is the detuning parameter. 
Another widely used model expands the non-linear term 
to give a cubic dependence on the field (third order Lamb 
theory ||), but this is not necessary here. Eq. (11) is 
written in a reference frame in which the frequency of the 
on steady state is zero 12 1. £(t) is a complex Langevin 
source term accounting for the stochastic nature of spon- 
taneous emission. It is taken as a Gaussian white noise 
of zero mean and correlations: 



(C(t)C(t'))=iA6(t-t') 



(12) 



where A measures the strength of the noise. 

By writing the complex variable E as E = x± + ix 2 and 
introducing a new dimensionless time such that t — » nt, 
the evolution equations become: 



X2 



(I 



b + x\ + x\ 
a 

b + x\ + x\ 



1 (xi -ax 2 )+^(t) (13) 



1 (a Xl + x 2 ) + &(*) (14) 



Where a = T/(k(3) and b = 1//3. £t(i) and are 
white noise terms with zero mean and correlations given 
by equation (||) with e = A/k. 

In the deterministic case (e = 0), these dynamical 
equations constitute a potential flow of the form (Q) 
where the potential V(x) is 

V(x 1 ,x 2 ) = - \ x i+ x l - a ln(b + xf + xj)] (15) 

and the matrix -D(x) (split into symmetric and antisym- 
metric parts) is: 



D=S+A= 



1 
1 



-a 
a 



(16) 



A simpler expression for the potential is given in Q and 
Jl7t valid for the case in which the gain term is expanded 
in Taylor series. 

According to our discussion above, the fixed points of 
the deterministic dynamics are the extrema of the poten- 
tial V(x): for a > b there is a maximum at (xi, x 2 ) = 
(corresponding to the laser in the off state) and a line of 
minima given by x\ + x\ = a — b (see Fig. 1). The asymp- 
totic stable situation, then, is that the laser switches to 



the on state with an intensity I = \E\ 



a — b. 



For a < b the only stable fixed point is the off state 1 = 0. 

In the transient dynamics, the symmetric matrix S is 
responsible for driving the system towards the line of 
minima of V following the lines of maximum slope of 
V. The antisymmetric part A (which is proportional to 
a) induces a movement orthogonal to the direction of 
maximum variation of V(x). The combined effects of 
S and A produce a spiraling trajectory in the (xi,x 2 ) 
plane. The angular velocity of this spiral movement is 
proportional to a. Asymptotically, the system tends to 
one of the minima in the line I = a — b, the exact loca- 
tion depending on the initial conditions. The potential 
decreases in time until it arrives at its minimum value: 
V(x\ + x\ = a - b) = -|(a ln(a) - a + 6). 

In the presence of moderate levels of noise, e ^ 0, the 
qualitative features of the transient dynamics remain the 
same as in the deterministic case. The most important 
differences appear near the stationary situation. As the 
final value of the intensity is approached and for a/0, 
the phase rotation slows down and the mean value of the 
phase 4> of the electric field E changes linearly with time 
also in the steady state, see Fig. 2. For a — there 
is only the ordinary phase diffusion around the circum- 



ference xf 



x\ 



b that represents the set of all 



possible deterministic equilibrium states [12[. Therefore, 
for a ^ the real and imaginary parts of E oscillate 
not only in the transient dynamics but also in the steady 
state, and while the frequency of the oscillations still de- 
pends on a (as well as e) , their amplitude depends on the 
noise strength e. 

We can understand these aforementioned features of 
the noisy dynamics using the deterministic Lyapunov po- 
tential V(xi, x 2 ). Since conditions 1-3 above are satis- 
fied, the stationary probability distribution is given by 
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( |io| ) with V(xi,X2) given by (|15|). By changing vari- 
ables to intensity and phase, we find that the probability 
density functions for / and are independent functions, 
P s t(4>) = 1/(2tt) is a constant and, 



Pst(I) 



-I/(2e) 



(b + I) 



a/(2e) 



(17) 



where the normalization constant is Z = 
(2e)* +1 e£r + 1, ^) and T(x,y) is the incomplete 
Gamma function. From this expression, we see that, 
independently of the value for e, P s t{I) has its max- 
ima at the deterministic stationary value I m — a — b. 
Starting from a given initial condition corresponding, for 
instance, to the laser in the off state, the intensity fluc- 
tuates around a mean value that increases monotonically 
with time. In the stationary state, the intensity fluctu- 
ates around the deterministic value I m = a — b but, since 
the distribution ( |l7| ) is not symmetric around I m , the 
mean value_(/) s t is larger than the deterministic value. 
By using ([L7]) one can easily find that 



(I) st = (a-b) + 2e 



exp(-b/2e){b/2e)z 

r(£ + i, 



£+1 



(18) 



An expression for the mean value of the intensity in the 
steady state was also given in [|l7| in the simpler case 
of an expansion of the saturation-term parameter in the 
dynamical equations. 

As mentioned before, in the steady state of the stochas- 
tic dynamics, the phase <j) of the electric field fluctuates 
around a mean value that changes linearly with time. Of 
course, since any value of <fi can be mapped into the in- 
terval [0, 27r), this is not inconsistent with the fact that 
the stationary distribution for <f> is a uniform one. We 
can easily understand the origin of this noise sustained 
Bow [po| : the rotation inducing terms, those proportional 
to a in the equations of motion, are zero at the line of 
minima of the potential V and, hence, do not act in the 
steady deterministic state. Fluctuations allow the sys- 
tem to explore regions of the configuration space (x%, X2) 
where the potential is not at its minimum value. Since, 
according to Eq. (plf), the mean value of / is not at the 
minimum of the potential, there is, on average, a non-zero 
contribution of the rotation terms producing the phase 
drift observed. 

The rotation speed can be calculated by writing the 
evolution equation for the phase of the electric field as: 



i+7- 1 ) B+ 77«" 



(W) 



where £(t) is a white noise term with zero mean value 
and correlations given by (||). By taking the average 
value and using the rules of the Ito calculus, one arrives 
at: 



(</>)= a 



b + I 



1 



(20) 



and, by using the distribution 
stochastic frequency shift: 



exp(-6/2e)(6/2e)^ 

1 ; T\ 

^ \2e ' x ' 2e) 



171), one obtains the 



(21) 



Notice that this average rotation speed is zero in the case 
of no detuning (a = 0) or for the deterministic dynamics 
(e = 0) and that, due to the minus sign, the rotation 
speed is opposite to that of the deterministic transient 
dynamics when starting from the off state. These results 
are in excellent agreement with numerical simulations of 
the rate equations in the presence of noise (see Fig. 2). 



III. CLASS B LASERS 

The dynamics of a typical class B laser, for instance 
a single mode semiconductor laser, can be described in 
terms of two evolution equations, one for the slowly- 
varying complex amplitude E of the electric field inside 
the laser cavity and the other for the carriers number 
N (or electron- hole pairs) JlCj ]. These equations include 
noise terms accounting for the stochastic nature of spon- 
taneous emission and random non-radiative carrier re- 
combination due to thermal fluctuations. Both noise 
sources are usually assumed to be white Gaussian noise. 

The equation for the electric field can be written in 
terms of the optical intensity / and the phase </> by defin- 
ing E = \/l e 1 ^ '. For simplicity, we neglect the explicit 
random fluctuations terms and retain, as usual fiot] , the 
mean power of the spontaneous emission. The equations 
are: 



dl 



(G(N, I) - 7)/ + A(3N 



i = 2( G ^ 



j)a 



^L=C- 7e N-G(N,I)I 
at 

G(N, I) is the material gain given by: 

g(N - N ) 



G(N,I) 



l + sl 



(22) 



(23) 



(24) 



(25) 



The definitions and typical values of the parameters for 
semiconductor lasers are given in Table 1. The first 
term of Eq. (^||) accounts for the stimulated emission 
while the second accounts for the mean value of the 
spontaneous emission power. Eqs. (^2] - |24|) are writ- 
ten in the reference frame in which the frequency of the 
on state is zero when spontaneous emission noise is ne- 
glected. The threshold condition is obtained by setting 
G(N, I) = 7, / = and neglecting spontaneous emission, 
i.e. = ^V + -- The threshold carrier injected per unit 
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time to turn the laser on is given by Cth = JeNth- Eq. 
( f23| ) shows that <j) 1S linear with N and slightly (due to 
the smallness of the saturation parameter s, see Table 1) 
nonlinear with I. 

Since in the deterministic case considered henceforth, 
the evolution equations for I and iV do not depend on 
the phase <fi, we can concentrate only in the evolution 
of / and N. One can obtain a set of simpler dimen- 
sionless equations by performing the following change of 
variables: 

V = —I, z = ^-(N-N ), T =lt (26) 
7 7 2 

The equations become then: 



^ = 2 

dr \ 1 + sy 

dz 



a-bz 
or 1 + sy 



1 | y + cz + d 

z y 



(27) 
(28) 

where we have defined a = -%{C — j e N ), b = 22=, 

c = d = 16 ^2 jV ° and s — j^. These equations form 
the basis of our subsequent analysis. The steady states 
are obtained by setting (p7|) and (28) equal to zero, i.e.: 



1 



//.,; : — — [ 2 (« - b) + d{l + bs) + cas + (29) 

(30) 



4(1 + bs) 1 
a(l + sy st ) 



b + y st (l + bs) 

where the constant v is given by: 

v = 4(<z - b) 2 + Ad(a + 6)(1 + bs) + d 2 {l + bs) 2 
+ c(8a + Aas(a + b) + 2das(l + bs)) + c 2 a 2 s 2 (31) 

There is another steady state solution for y st given by Eq. 
( p9| ) (with a minus sign in front of y/v) which, however, 
does not correspond to any possible physical situation, 
since y s t < 0. For a value of the injected carriers per unit 
time below threshold (C < Cth, equivalent to a — b < 0) 
y s t is very small. This corresponds to the off solution 
in which the only emitted light corresponds to the spon- 
taneous emission. Above threshold, stimulated emission 
occurs and the laser operates in the on state with large 
y s t- In what follows, we will concentrate in the evolution 
following the laser switch-on to the on state. 

It is known that the dynamical evolution of y and z 
is such that they both reach the steady state by per- 
forming damped oscillations [0 whose period decreases 
with time. This fact is different from the usual relaxation 
oscillations that are calculated near the steady state by 
linearizing the dynamical equations. The time evolution 
of y and z is shown in Fig. 3a for some parameters (for 
another values of the parameters equivalent results are 
obtained), while the corresponding projection in the y, 
z phase-plane is shown in Fig. 4. We are interested in 
obtaining a Lyapunov potential that can helps to explain 



the observed dynamics. This study was done in (Tj] with- 
out considering neither the saturation term, nor the mean 
value of the spontaneously emission power, and under 
those conditions an expression for the period of the tran- 
sient oscillations was obtained. In our work, we calculate 
the period of the oscillations by taking into account these 
two effects. The period is obtained in terms of the po- 
tential, by assuming that the latter has a constant value 
during one period. It will be shown that this assumption 
works reasonably well and gives a good agreement with 
numerical calculations. Near the steady state, the relax- 
ation oscillations can also be calculated in this form, but 
the potential is almost constant and consequently so is 
the period. 

The evolutions equations (p7f), ( p8|) can be casted in 
the form of Eq. (^) with the following Lyapunov poten- 
tial: 

Aa 1 

V(y, z) = Ax y + A 2 y 2 + A 3 hx(y) + — + - B 2 (y, z) 

V 2 

(32) 



where 



A-i = as + bs sd (1 + bs) as 2 c 

2 2 4 4 

A 2 = S - (1 + bs) 

A 3 = ~ [a - b + (ac + bd) s + |] 



A 4 



(ac + bd) 



„/ , (d + cz) , 

B{y, z) = z-l-sy+ i— '- 1 + sy) . 

2y 

The corresponding (non-constant) matrix D is given by: 



D 





d\2 



-di2 
d 2 2 



being 



dl2 



dii = 



Ay 2 



(1 + Sy) [2y + c (1 + Sy)] 
Ay [(1 + 2s + bs) y 2 + by - 



d + cz] 



(l + sy) [2y + c{l + sy)Y 



(33) 

(34) 
(35) 



This potential reduces to the one obtained in ref. JTJ] 
when setting c — d — s — (which corresponds to set the 
laser parameters (5 = s = {)). As expected, non- vanishing 
values for the parameters s and (3 increase the dissipative 
part of the potential (^22), associated with the damping 
term. This result was pointed out in [glj when linearizing 
the rate equations around the steady state. 

The equipotential lines of (32) are also plotted in Fig. 
4. It is observed that there is only one minimum for V 
and hence the only stable solution (for this range of pa- 
rameters) is that the laser switches to the on state and 
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relaxes to the minimum of V. The movement towards 
the minimum of V has two components: a conservative 
one that produces closed equipotential trajectories and a 
damping that decreases the value of the potential. The 
combined effects drives the system to the minimum fol- 
lowing a spiral movement, best observed in Fig. 4. 

The time evolution of the potential is also plotted in 
Fig. 3b. In this figure it can be seen that the Lyapunov 
potential is approximately constant between two consec- 
utive peaks of the relaxation oscillations (This fact can 
be also observed with the equipotential lines of Fig. 4). 
This fact allows us to estimate the relaxation oscillation 
period by approximating V(y, z) = V, constant, during 
this time interval. When the potential is considered to 
be constant, the period can be evaluated by the stan- 
dard method of elementary Mechanics: z is replaced by 
its expression obtained from ( p7|) in terms of y and y 
(the dot stands for the time derivative) in V(y, z). Using 
the condition that V(y, z) — V — constant, we obtain 
an equation for y of the form: F (y, y) = V. From this 
equation, we can calculate the relaxation oscillation pe- 
riod (T) by integrating over a cycle. This leads to the 
expression: 



where: 



T 



1 



sy 



dy 



-1MV2 



,„ V [2{V - A iy - A 2 y 2 - A 3 hx(y) - A iy -^)} 

(36) 



where yo and y\ are the values of y that cancel the de- 
nominator. We stress the fact that the only one approx- 
imation used in the derivation of this expression is that 
the Lyapunov potential is constant during two maxima 
of the intensity oscillations. The previous equation for 
the period reduces, in the case c = d = s = 0, to the one 
previously obtained by using the relation between the 
laser dynamics and the Toda oscillator derived in [pr[ . 
Evaluation of the above integral shows that the period 
T decreases as the potential V decreases. Since the Lya- 
punov potential decreases with time, this explains the 
fact that the period of the oscillations in the transient 
regime decreases with time. In Fig. 5 we compare the 
results obtained with the above expression for the period 
with the one obtained from numerical simulations of the 
rate equations (j27l), (p8|). In the simulations we compute 
the period as the time between two peaks in the evolution 
of the variable y. As can be seen in this hgure, the above 
expression for the period, when using the numerical value 
of the potential V, accurately reproduces the simulation 
results although it is systematically lower than the nu- 
merical result. The discrepancy is less than one percent 
over the whole range of times. 

It is possible to quantify the difference between the ap- 
proximate expression ( p6| ) and the exact values near the 
stationary state. In this case expression (R6I) reduces to : 



2 7T 



VEF - H 2 



(37) 



1 A a 



Aa 



E = 2 A 2 --^ + ^ + - 



2y, 



F 
H 



1 + c 



(1 + sy st ) 



Vst 
n 2 



(d + cz st ) 



2y s t 
c(l + sy st ) 



2y s 



s + 



(d + cz st ) 

2y 2 st 



and d±2,st is the coefficient d\ 2 calculated in the steady 
state. The period of the relaxation oscillations near the 
steady state can be obtained by linearizing eqs. ( |27j ) and 
( p8| ) after a small perturbation is applied. The frequency 
of the oscillations in the steady state is the imaginary 
part of the eigenvalues of the linearized equations. This 
yields a period: 



2 7T 



di 2 , flt vef - m 



d 2 

(I 90 



-,-1/2 



d\ 2 ,st KEF 



H 2 ) 



(38) 



The difference between (J37]) and ( |38| ) vanishes with 
d 2 2,st (i-e. d 22 in the stationary state). Since EF — H 2 is 
always a positive quantity, our approximation will give, 
at least asymptotically, a smaller value for the period. 

In order to have a complete understanding of the vari- 
ation of the period with time, we need to compute the 
time variation of the potential V(r) between two con- 
secutive intensity peaks. This variation is induced by 
the dissipative terms in the equations of motion. We 
have not been able to derive an expression for the vari- 
ation of the potential (see [ fL4| for an approximate ex- 
pression in a simpler case). However, we have found that 
a semi-empirical argument can yield a very simple law 
which is well reproduced by the simulations. We start 
by studying the decay to the stationary state in the lin- 
earized equations. By expanding around the steady state: 
y = Vst + Sy, z = z st + 5z, the dynamical equations imply 
that the variables decay to the steady state as: Sy(r), 
6z(t) oc exp(— |r), where: 



p = d 



22, st 



F 



(39) 



Expanding V(y, z) around the steady state and taking 
an initial condition at tq we find an expression for the 
decay of the potential: 



In [V(t) - V st ] = In [V(t ) - V st ] - p (i 



(40) 



In Fig. 6 we plot ln[V (r)— V s t] versus time and compare it 
with the approximation ([hj]) . One can see that the latter 
fits ln[y(r) — Vst] n °t only near the steady state (where 
it was derived), but also during the transient dynamics. 
The value of To, being a free parameter, was chosen at 
the time at which the first peak of the intensity appears. 
Although other values of tq might produce a better fit, 



G 



the one chosen here has the advantage that it can be 
calculated analytically by following the technique of ref. 
p2| . It can be derived from Eq. (|6|) that the period T is 
linearly related to the potential V . This, combined with 
the result of Eq. ([lO]), suggests the semi-empirical law 
for the evolution of the period of the form: 

In [T(r) - T st ] = In [T(r ) - T st ] - p (r - r ) (41) 

This simple expression fits well the calculated period 
not only near the steady state, but also in the transient 
regime, see Figs. 5 and 7. The tiny differences observed 
near the steady state are due to the fact that the semi- 
empirical law, Eq. (|4l]) , is based on the validity of rela- 
tion Eq. ( p6|) between the period and the potential. As 
it was already discussed above, that expansion slightly 
underestimates the asymptotic (stationary) value of the 
period. By complementing this study with the procedure 
given in p3] to describe the switch-on process of a laser, 
and valid until the first intensity peak is reached, we can 
obtain a complete description of the variation of the os- 
cillations period in the dynamical evolution following the 
laser switch-on. 



IV. SUMMARY 

In this work we have used Lyapunov potentials in the 
context of laser dynamics. For class A lasers, we have 
explained qualitatively the observed features of the de- 
terministic dynamics by the movement on the potential 
landscape. We have identified the relaxational and con- 
servative terms in the dynamical equations of motion. In 
the stochastic dynamics (when additive noise is added 
to the equations), we have explained the presence of a 
"noise sustained flow" for the phase of the electric field 
as the interaction of the conservative terms with the noise 
terms. An analytical expression allows the calculation of 
the phase drift. 

In the case of class B lasers, we have obtained a Lya- 
punov potential valid only in the deterministic case, when 
noise fluctuations are neglected. We have found that the 
dynamics is non-relaxational with a non-constant matrix 
D. The fixed point corresponding to the laser in the on 
state is interpreted as a minimum in the potential land- 
scape. By observing that the potential is nearly constant 
between two consecutive intensity peaks during the re- 
laxation process towards the steady state, but still in 
a highly non-linear regime, we were able to obtain an 
approximate expression for the period of the oscillations. 
Moreover, we have derived a simple exponential approach 
of the period of the oscillations with time towards the pe- 
riod of the relaxation oscillations near the steady state. 
This dependence appears to be valid after the first inten- 
sity peak following the switch-on of the laser. A possible 
extension of our work could be to consider the presence 
of an external field, which is numerically studied in (B3| . 
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Table 1 



PARAMETERS 




VALUES 


C 


Carriers injected per unit time. 


> threshold 


1 


Cavity decay rate. 


0.5 ps~ L 




Carrier decay rate. 


0.001 ps~ l 


No 


Number of carriers at transparency. 


1.5 x 10 s 


9 


Differential gain parameter. 


1.5 x 10~ 8 ps- L 


s 


Saturation parameter. 


io- s - io-'' 




Spontaneous emission rate. 


io^v -1 


a 


Lincwidth enhancement factor. 


3-6 



8 



FIGURE CAPTIONS 



Fig. 1: Potential for a class A laser, Eq. ( |l5| ) with the parameters: a = 2, b = 1. Dimensionless units. 



Fig. 2: Time evolution of the mean value of the phase <f> in a class A laser, in the case a = 2, 6 = 1, e = 0.1. For a = 
(dashed line) there is only phase diffusion and the average value is for all times. When a = 5 (solid line) there is a 
linear variation of the mean value of the phase at late times. Error bars are incluted for some values. The dot-dashed 
line has the slope given by the theoretical prediction Eq. (^l|). The initial condition is taken as x\ — X2 — and the 
results were averaged over 10000 trajectories with different realizations of the noise. Dimensionless units. 



Fig. 3: a) Normalized intensity, y (solid line) and normalized carriers number, z/40 (dot-dashed line) versus time 
in a class B laser obtained by numerical solution of Eqs. ( p7| ) and fl28|). b) Plot of the potential ([52|). Parameters: 
a = 0.009, b = 0.004, s = 0.5, c = 3.2 x 10~ 9 , d = 1.44 x 10~ 8 which correspond to physical parameters in Table 1 
with C = 1.2Cth- The initial conditions are taken as y = 5 x 10~ 8 and z = 0.993. Dimensionless units. 



Fig. 4: Number of carriers versus intensity (scaled variables). The vector field and contour plot (thick lines) are also 
represented. Same parameters than in Fig. 3. Dimensionless units. 



Fig. 5: Period versus time in a class B laser. Solid line has been calculated as the distance between two peaks of 
intensity, with triangles plotted at the begining of each period; dashed line has been calculated using the expression 
(|36|), with the value of the potential V obtained also from the simulation; dotted line corresponds to the semi-empirical 
expression (|4l|). Same parameters than in Fig. 3. We have used tq = 55.55, coinciding with the position of the first 
intensity peak. Dimensionless units. 



Fig. 6: Logarithm of the potential difference versus time in a class B laser (solid line), compared with the theoretical 
expression in the steady state ( [40| ) (dashed line). Same parameters than in Fig. 3 and t as in Fig. 5. Dimensionless 
units. 



Fig. 7: Logarithm of the period difference versus time in a class B laser. Triangles correspond to the period calculated 
from the simulations as the distance between two consecutive intensity peaks, at the same position than in Fig. 6. 
The dashed line is the semiempirical expression Eq. (p|). Same parameters than in Fig. 3 and To as in Fig. 5. 
Dimensionless units. 
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Mayol et al. Fig. 1 




Mayol et al. Fig. 2 
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Mayol et al. Fig. 4 
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Mayol et al. Fig. 5 




Mayol et al. Fig. 6 
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Mayol et al. Fig. 7 



13 



